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Abstract 

■ The interacting dark energy model with interaction term Q = 

\ m H p m + XdH Pd is considered. By studying the model near the transi- 

• tion time, in which the system crosses the u> = — 1 phantom-divide-linc, 

| the conditions needed to overcome the coincidence problem is inves- 

ts i tigated. The phantom model, as a candidate for dark energy, is con- 

sidered, and for two specific examples, the quadratic and exponential 
phantom potentials, it is shown that it is possible the system crosses 
the ijj = — 1 line, meanwhile the coincidence problem is alleviated, the 
two facts that have root in observations. 

cr 
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1 Introduction 



Nowadays based on astrophysical data it is believed that the universe is 
accelerating [1]. The origin of this acceleration is still unknown and differ- 
ent models have been proposed to elucidate this subject. One picture is 
the assumption that nearly 70% of the universe is composed of a smooth 
energy component with negative pressure dubbed as dark energy. A sim- 
ple candidate for dark energy is the cosmological constant [2] which suffers 
from conceptual problems such as fine-tuning and coincidence problems [3]. 
Therefore alternative models, e.g., introducing dynamical exotic fields such 
as scalar fields with suitably chosen potentials, have been introduced [4]. 

In dark energy models, the ratio of matter to dark energy density, r, is 
expected to decrease rapidly (proportional to the scale factor) as the uni- 
verse expands, but observations show that these densities are of the same 
order today. To solve this problem (known as coincidence problem), one can 
adopt an evolving dark energy field with suitable non-gravitational interac- 
tion with matter [5], [6]. Various models corresponding to different forms of 
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interaction, leading to a constant or slowly varying (soft coincidence) r at 
late times, have been proposed [7]. 

Some present data seems to favor an evolving dark energy, corresponding 
to an equation of state (EOS) parameter less than uj = — 1 at present epoch 
(phantom regime) from u > — 1 in the near past (quintessence regime) [8]. 
So another cosmological coincidence problem may be proposed: why to = — 1 
crossing is occurred at the present time [9] . 

In [10], it was shown that uj = — 1 crossing in models including matter 
and phantom scalar field is either impossible or unstable with respect to 
cosmological perturbations. However, this transition may be possible for 
scalar-tensor theories [11], multi-field models [12], and coupled dark energy 
models with specific couplings [13], [14]. 

In [15], the transition from quintessence to the phantom phase in the 
quintom model was considered in the slow roll approximation. By studying 
the Friedman equations near the transition time, it was shown that in non- 
interacting quintom model, r ~ at transition time. This lies in the fact 
that the main part of the dark energy at transition time corresponds to the 
quintom potential. By considering interaction between cold dark matter and 
dark energy, the mutual energy exchange between two fluids will be allowed 
and the coincidence problem may be alleviated. 

In this paper we consider dark energy model composed of a phantom 
scalar field interacting with cold dark matter. We try to elucidate the con- 
nection between the coincidence problem and uj = — 1 crossing (second cos- 
mological coincidence problem). 

It may be worth noting that the phantom models suffer from the quan- 
tum instability problem. Because the phantom fields have negative kinetic 
energy, it is possible that a phantom particle decays into an arbitrary number 
of phantoms and ordinary particles, such as gravitons. It can be shown that 
the decay rates of these interactions are infinite, which indicates that the 
phantom models are dramatically unstable. But if we think of these models 
as the low-energy effective theories, with the fundamental fields having pos- 
itive kinetic energy, then we should use a momentum cutoff A in calculating 
the decay rates. In this way, it can be shown that, for A ~ M p \, the lifetimes 
can become larger than the age of the universe when one chooses suitable 
phantom-gravity interaction potentials, and this removes the quantum in- 
stability of these kinds of phantom models [16]. 

The Scheme of the paper is as follows. After the Introduction, we con- 
sider the dark energy model with interaction term Q = \ m Hp m + XjHpd 
in section two. By restricting ourselves to times t « 1 around the tran- 
sition time ( ho is the Hubble parameter), we study the general properties 
of interacting dark energy models and the necessary conditions needed to 
cross uj = — 1 line are obtained. These results are insensitive to the origin of 
the dark energy. In section three we assume that dark energy is composed 
of phantom scalar field interacting with cold dark matter. After a general 
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discussion, we illustrate, via two specific examples, how the necessary con- 
ditions for uj = — 1 crossing can alleviate the coincidence problem. It is seen 
that it is possible to tune the parameters such that r = 3/7 at transition 
time. 

We use units fi = c = G = 1 throughout the paper. 



2 uj = — 1 crossing in interacting dark energy model 

We consider a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe containing dark energy and dark matter fluids. In terms of dark 
energy density pd and matter energy density p m , the Hubble parameter is 
given by Friedmann equation 

H 2 = —p = —{p m + Pd), (1) 

where p is the total energy density. By introducing = pd/p and Q m = 
Pm/p, eq.(l) can be written as Qd+f2 m = 1, which indicates that the universe 
is spatially flat. The derivative of the Hubble parameter with respect to the 
comoving time can be extracted from Einstein equations. The result is 

H = -4ir(p d + P d + Pm ). (2) 

Pd is the pressure of the dark energy fluid and the dark matter is assumed 
to be pressureless. The equation of state of the universe is P = up, where 
P = Pd is the pressure and uj is the equation of state parameter which can 
be written as 

For an accelerated universe we have uj < —1/3. When —1 < uj < —1/3, 
the universe is in quintessence phase and when uj < — 1, the universe is in 
phantom phase. In the following, we assume that the dark matter and dark 
energy components can interact through the following source term: 

Q = ^mHpm + XdHpd, (4) 

where X m and are two real constants. For special choices such as X m = 0, 
Xd = or Xd = X m , eq.(4) reduces to the interaction terms which have been 
considered before [5]. The other forms of interaction terms, not necessarily 
suitable for our purpose, have been also considered in the literature [17]. 

Because of the interaction term, we have not the conservation of partial 
stress-energy tensors of matter and dark-energy: T,^K , v = —T^. v / 0. 
In fact, the projection of this nonconservation equation along the velocity 
of the whole (comoving) fluid U u (which was taken to be the same as the 
velocities of each of the fluid components) is [6] 

U v Tf\ = -U V T,% =-Q. (5) 
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Note that the coupling (4) can be written as a scalar as follows 

Q = \u,U v {\ m T { ^ + \ d T$)U« a . (6) 

For FLRW metric, the equation (5) reduces to: 

p d + 3H{p d + P d ) = -Q, 

p m + 3Hp m = Q. (7) 

Using eq.(l), eq.(7) can be written as 

3 o 

p d + (3 + X d - X m )Hp d + 3HP d = -—X m H 3 , 

Pm + (3 + X d - X m )Hp m = ^A d # 3 . (8) 

Using eq.(8), the evolution equation of the ratio of energy densities of dark 
matter and dark energy, denoted by r = p m /p d , reads 

r = r(r + 1) Uuj + X m + ^\ H. (9) 



From 



eq.(9) then results in 



n d = — , (io) 



- 1 Xg^d X rn 

u ~ 3H i - n d 3(1 - n d ) 3 ' ( ' 

In the vicinity of transition time from quintessence to phantom era, 
U > —1 goes to uj < —1, so H must change sign from H < to H > 0. At 
transition time we have H = and lo = — 1. The dark energy equation of 
state parameter uj d is defined through P d = io d Pd- Therefore ujp = u d p d or 
Q d ^d = u. Using 

3u: d H 2 + 2H + 3H 2 
Pm = — 



8truj d 
2H + 3H 2 

Pd = 5 , (12) 

and eq.(8), one can obtain the following equation for the Hubble expansion 
H + (6 + X d - X m + 3u d )HH + - [(3 - X m )oj d + 3 + X d - X m ] H 3 
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For a constant Ud, we arrive at the result of [18]. At H = we obtain 



H = ~[(Z-\ m )u d + ?> + \ d -\ m ]H\ (14) 

Note that H > 0, therefore for a constant u>d, the sign of H does not change. 
This shows that in constant-u;^ approximation, the system can cross the 
u = — 1 line only once. This is because the transition from quintessence to 
phantom phase needs positive H (at transition time), while the vice versa 
needs negative H. In the following, we consider u>d as a function of time. 
At transition time, we obtain from eq.(ll) 

^ = (3 . Xm+Xd) ^ + _^_ y (15) 



which results in 



Insertion of 



(3 - A m + X d )n d + ^ > 0. (16) 



n d = 



8?r p d 
3 # 2 ' 



into eq.(16), leads to 



2H 

(3 - A m + Adjffpd + p d - -j^Pd > 0, (18) 

at transition time. 

We assume that, in the neighborhood of transition time, the Hubble 
parameter is a differentiable function of time. The Taylor expansion of H 
at transition time, which we take as t = 0, can be written as [15] 

H = ho + htt + 0(t a+1 ), a>2,hx (19) 

foo = H(t = 0), a is the order of the first nonzero derivative of the Hubble 
parameter at transition time, and hi = (1 / 'a\)d a H / 'dt a \t=Q. The transition 
from quintessence to phantom phase occurs if a is an even positive integer 
and hi > 0. We also consider the following expansions for p m , and pd 
at t = 

V d = uo + uit? + 0(t l3+1 ), 

Pm = PmO + Pmlf 1 + 0(t 1+1 ), 

Pd = Pdo + Pdit + 0(t e+1 ), (20) 
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respectively. (3, 7, and 9 are the orders of the first nonzero derivatives of 
Qd, Pmi an d pd at t = 0, respectively. Note that the above expansions are 
valid until t << h^ 1 , which is completely reasonable since 1 is of order of 
age of our universe. 

To obtain the relation between the parameters a, (3, 7, and 6, we proceed 
as follows. For (3 7^ 1, if we expand both sides of eq.(ll) at t = 0, the first 
resulting term of the right hand side, with nonvanishing power of t, is t 13 - 1 
while the left hand side (after t°) begins with t a ~ l . So if (3 7^ 1, we must 
have a = (3. In this case, (3/87r)H 2 Qd = pa results = (3(= a). For (3 = 1, 
this equation results in 8 = (3{= 1). Therefore always 6 = (3. From eq.(12) 
it is clear that, for a constant ujd, we must have (3 = a — 1 which leads to 
(3=1 and a = 2. 

In the case (3^1, comparing the coefficients of i°-terms of eq.(ll) gives 

uo — - — ^ m (21) 
Ad — A m + 3 

and equating the coefficients of t a ~ l -terms results in 

2(1 -n ) V ' 

This relation shows that the transition is possible only if u\ > 0. 
In the case (3=1, the same procedure leads to 

u\ = (A m - Ad - 3)h u + (3 - X m )h , (23) 

(Ad - A m + 3)/i n fe + {k + l)u k+1 = 0, 1< fc < a - 2, (24) 

and 

_ (Ad - A m + 3)/ign Q _i + a/i- n a 
1 ~ 2a(l - n ) ' 

The Taylor expansion of r at £ = is 

r = r + r 1 t (S + O(t 2 ), (26) 

where ro = nQ 1 — 1 and n = —ui/(uq). 



3 Interacting phantom dark energy model and co- 
incidence problem 

In this section we assume that the origin of the dark energy is a phantom 
scalar field <p. So 

Pd = ~^> 2 + V(<f>), 

Pd = -\P-V{4>), (27) 
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where V(4>) > is the phantom potential, ujd is given by 

Wd ~ -^ + vw (28) 

therefore oj d < — 1- For fl d uj d < — 1> the universe is in the phantom and for 
Q d uJd > —1 it is in the quintessence phase. 
The field equation of (f> is 

■ dV 

(/)(</> + 3H<f> - — ) = Q. (29) 

This can be derived by putting eq.(27) back into eq.(7). From eq.(27) we 
obtain 

<p 2 = -{l+u d )pd 
2V{4>) = (l-u d ) Pd . (30) 

The second equation of (30) can be written as 

• dV~\y) 



dy 



-y, (31) 



where y = (1 — ujd)p d /^-, and V 1 is the inverse function of V. Eq.(31) and 
the first equation of (30) then lead to 

which after some calculation can be rewritten as 

i 2 fdV~ 1 (v)\ 2 4 
Sl d -u-3H(l+u)(n d -u>) , {U) ) =—(n d + u). (33) 

1 \ dy J p 

This equation together with our previous results may be served to find some 
necessary conditions for uj = — 1 crossing in interacting phantom dark energy 
models, including the domain to which no belongs. We will try to obtain 
a relation between the coincidence problem and the behavior of the system 
at transition time. For example, for the case (3 = 1, if one obtains h\ as a 
polynomial of uq, then restricting h\ to positive values, which is necessary 
for transition, will restrict the value of uq to a subset of (0, 1). By choosing 
the appropriate parameters, then it becomes possible to prevent r to be or 
very large. For (3^1 cases, eq.(21) determines the value of r at transition 
time, which again can be chosen to be O(l). In this way the occurrence 
of a; = —1 crossing and the alleviation of the coincidence problem can be 
achieved simultaneously. 

In the following, we will show these points via some specific examples. 
In these examples we restrict ourselves to the case a = 2. 
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3.1 Phantom field with square power law potential 

For V{4>) = (l/2)m 2 cf) 2 , eq.(33) becomes 



2mJto 2 - n 2 d 



n d -LO-3H(l+uj)(n d -Lo) . (34) 



In the following, we adopt that, in the quintessence phase and near the 
transition time, Q d > or equivalently r\ < [19]. Therefore 



2m^Lo 2 -tt 2 = n d -u-3H(l+u)(n d -u}) (35) 

can be used in the neighborhood of transition time. Taking f3 = 1 (the case 
(3 ^ 1 will be discussed later), the expansion of eq.(35) at t = then results 
in 

2 " 1 ^ t + 0{t) 

+ H + ( 2 " 2+4 | +4 ^) 4 + ° (42) - (36) 



As a result we arrive at 

4/ii 

3 



The necessity of quintessence to phantom phase transition, i.e. /ii > 0, then 
results in 

«i > 2m^l- u 2 . (38) 
Using eq.(23), we can write the above inequality in terms of uq: 

au + b < cyjl- Uq. (39) 

We have defined a = X m — X d — 3, b = 3 — X m , and c = 2m/ho. 

To study the solutions of eq.(39), we consider two situations. The first 
possibility is: 

au + b < 0, (40) 

which leads to r\ > 0. This conflicts with the assumption r\ < 0, or 
equivalently fl d > at transition time, and therefore is not acceptable. 
The second possibility is auo + b > which leads to 

V(uq) := (a 2 + c 2 )ul + 2abu + b 2 - c 2 < 0. (41) 

If a 2 — b 2 + c 2 < 0, V has no real roots and its sign does not change. But 
V(l) > 0, therefore eq.(41) is not satisfied in this case. For a 2 — b 2 + c 2 > 0, 
V has two roots which we denote by ur\ and ur2- Eq.(41) is satisfied if the 
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value of fid at transition time is restricted to the intersection of the intervals 
(uri,u R2 ) and (0, 1) 

u e{0,l)f){u R1 ,u R2 ). (42) 

So if {uri,u R 2) C (0,1), by choosing the appropriate parameters a, b, and 
c, we can obtain the desired order of magnitude: ~ 0(1) for ro = l/«o — 1. 
The Sturm sequences at and 1 are: 



S(0) = 



a 2 + c 2 



and 



S(l) 



(a + b) 2 , 2(a 2 + c 2 + ab), 



c 2 (a 2 -b 2 + c 2 ) 
a 2 + c 2 



(43) 



(44) 



Using Sturm theorem, one can show that for 

a 2 - b 2 + c 2 > 0, 
a 2 + c 2 + ba> 0, 
6 2 - c 2 > 0, 
ab < 0, 



the two roots of V belong to (0, 1). In this way we have 



-ab — c\J a 2 + c 2 — b 2 —ab + c\/ a 2 + c 2 — b 2 

— < ild < — 



a 2 + c 2 



a 2 + c 2 



(45) 



(46) 



at transition time. As an example, consider the case X m = 1, = 2 and 
c = 1. In this case ur\ = 0.258 and u R 2 = 0.682, therefore 0.46 < r < 2.8 
at transition time. Note that (u R i,ur2) may be more tightened by choosing 
appropriate o, b and c, e.g. for c = 1, a = 7, and b = —5, which correspond 
to X m = 8 and = —2, we have uri = 0.6 and u R 2 = 0.8, therefore 
0.6 < uq < 0.8 in agreement with the expected value uq = 0.7 and ro = 3/7. 

In P ^ 1 cases, we have /3 = a. For a = 2, following the same method 
ending to eq.(36), we expand both sides of eq.(35). It is obtained, up to the 
order 0(t), 



m 



1 



2hi 



(47) 



which results in h\ > 0. uq and u\ are given by eqs.(21) and (22), respec- 
tively. Higher orders of t in the Taylor expansion of eq.(35) determine the 
other coefficients of the Taylor expansion of H an f^. Since h\ > induces 
no additional constraint on no, the appropriate parameters a, 6, and c can 
lead to the desired values for ro. For a > 2, the aforementioned expansion 
leads to m\J\ — Uq = which is ruled out by the assumption that uq ^ 1. 
Therefore in (3 7^ 1, the choice a = 2 is the only eligible one. 
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Eq.(35) with Cl^ < and eq.(34) with minus sign can be also investigated 
by the same method. In brief, it is shown that, in interacting phantom model 
with V(4>) = (l/2)m 2 (p 2 phantom potential and the interaction Q-term of 
eq.(4), it is possible to choose the parameters such that both the u> = —1 
crossing and r = 0(1) occur. In special case which leads to eq.(47), one 
can tune the parameters such that tq has no choice but the desired value 
3/7. 

3.2 Exponential potential 

Consider the following potential 

V = v exp(\(f>), A > 0, v > 0. (48) 
Eq.(33) then results in 

H\(io 2 -Q 2 d )^(n d -Lo)^ = ± \n d - u - 3H(1 + u){Q d - u)] , (49) 



where A = a/3/ (8tt)X. As the previous example, we consider the upper sign 
of eq.(49) which is a result of the assumptions > and u < in the 
vicinity of transition time. 

By Taylor expansion of the both sides of eq.(49) at transition time, we 
obtain the following equation for a = 2 and (5 = 1 

1J| = -auo-b + X(l-u 2 )^(l + u )^. (50) 

a and b are defined by the same relations as for the first example. h\ > 
then results in 

au + b< A(l-^)5(l + u )i (51) 

Eq.(51) is satisfied in two cases: (i) auo + b < 0, which makes r\ negative, 
or firf < 0, and is not acceptable, (ii) auQ + b > 0. In this case we must have 

ul + (A 2 + l)ul + {2AB - l)u + B 2 - 1< 0, (52) 

where A = a/(A 2 ) and B = b/(\ 2 ). We also assume B 2 > 1. Eq.(52) is 
satisfied only if the polynomial 

Q(u ) := u% + (A 2 + l)ul + (2AB - l)u + B 2 - 1 (53) 

has real roots. Following Descartes rule, B 2 — 1 > and 2AB — 1 < are 
necessary conditions for Q(uq) to have two real positive roots. The domain 
to which uo in eq.(52) belongs is the intersection of (0,1) and (uri,ur2), 
where uri and ur2 are the roots of Q(uq). So by appropriate choosing 
of A and B, one can restrict uq to the domain allowed by astrophysical 
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data. To do so, we construct the Sturm sequence corresponding to the cubic 
polynomial Q(uq) at and 1. They are 



S(0) 



and 



S(l) = 



1, 2AB - 1, -(2AB - l)(A 2 + 1) + (1 - B 2 
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B' 
9 

4 (A 4 + 2A 2 - 6AB + 4) 2 



D 



(54) 



A 2 + 2AB + B 2 , 2 A 2 + 2AB + 4, -(2A A + 2A 3 5 + 3A 2 
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-10 AB - 9B' + 16), - 



D 



4 (A 4 + 2A 2 - 6,45 + 4) 2 



(55) 



In above equations, D > is discriminant of the polynomial Q(tio)- By 
implying Sturm theorem, it can be verified that, in order to have two real 
roots in interval (0, 1), the parameters A and B must satisfy, besides the 
previous mentioned conditions B 2 — 1 > and 2AB — 1 < 0, the following 
inequalities 



A + AB + 2 > 0, 

2A 4 + 2A 3 B + 3A 2 - 10AB - 9B 2 + 16 > 0. 



(56) 



For example, for A = 1, X m = 1, and A^ = 2, we obtain 0.23 < uq < 0.73 
which is in agreement with uq ~ 0.7 obtained from astrophysical data. 

Now we consider (5^1. For a = 2, the Taylor expansion of the both 
sides of eq.(49), with upper sign, result in 

\ -^2 - A v 7 ! - -"oVl + u + 
3 ho 




4 T ^-2n 2 + 4 



^{1 + uq) 2 A/iiVl - V 4 A/nvT+tt 



/i 



^ = A(1 



3 /igyrT^ 3 ^yr 




;(i + u ) 



(58) 



which implies h\ > 0. By suitable choosing of A m and A^, one can obtain 
the appropriate value for r at transition time. For a > 2, the aforemen- 
tioned expansion leads to A(l — Uq) 5 (1 + no) 5 = which is ruled out by the 
assumption that uq / 1. Therefore a = 2 is the only allowed case for (5 / 1. 
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4 Conclusion 



In this paper, by considering the energy exchange between cold dark mat- 
ter and dark energy (see eq.(4)), we study the possibility of simultaneous 
occurrence of two phenomena, the coincidence problem, and uj = — 1 cross- 
ing, from w > -1 to w < -1. We consider the physical quantities near 
the transition time t = 0, through eqs.(19) and (20). The transition occurs 
for positive h± and even a, the parameters which has been introduced in 
eq.(19). 

The equation of state parameter ui is expressed by eq.(ll) and the po- 
tential of phantom field, as a candidate of dark energy, enters in eq.(33). 
We studied the perturbative solutions of these equations, near t = 0, for two 
specific potentials, i.e. the quadratic and exponential potential. It is shown 
that always 6 = (3, and for (3 ^ 1, a = f3 (see eq.(20) and its subsequent 
discussion). For a = 2, as a first acceptable solution for uj = — 1 crossing, 
it is shown that, in both examples, it is possible to choose the parameters 
such that, besides the satisfaction of dynamical equations, the occurrence 
of a; > -1 to w < —1 transition allows the ratio ro = (p m / Pd)t=o to be 
around the desired value 3/7. This proves the possibility of solving these 
two problems in a unique framework. 

For a > 2 and [3 = 1, it can also be shown that it is possible to choose 
the parameters such that the above mentioned properties are achieved. 

Acknowledgement: We would like to thank M. Honardoost for useful 
discussion. This work was partially supported by the "center of excellence 
in structure of matter" of the Department of Physics. 
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